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Abstract. 

Testing deviation of GR is one of the main goals of the proposed Laser Interferometer Space 
Antenna. For the first time, we consistently compute the generation of gravitational waves 
from extreme-mass ratio inspirals (stellar compact objects into supermassive black holes) in a 
well-motivated alternative theory of gravity, that to date remains weakly constrained by double 
binary pulsar observations. The theory we concentrate on is Chern-Simons (CS) modified 
gravity, a 4-D, efi'ective theory that is motivated both from string theory and loop-quantum 
gravity, and which enhances the Einstein-Hilbert action through the addition of a dynamical 
scalar field and the parity-violating Pontryagin density. We show that although point particles 
continue to follow geodesies in the modified theory, the background about which they inspiral 
is a modification to the Kerr metric, which imprints a CS correction on the gravitational waves 
emitted. CS modified gravitational waves are sufficiently different from the General Relativistic 
expectation that they lead to significant dephasing after 3 weeks of evolution, but such dephasing 
will probably not prevent detection of these signals, but instead lead to a systematic error in 
the determination of parameters. We end with a study of radiation-reaction in the modified 
theory and show that, to leading-order, energy-momentum emission is not CS modified, except 
possibly for the subdominant effect of scalar-field emission. The inclusion of radiation-reaction 
will allow for tests of CS modified gravity with space-borne detectors that might be two orders 
of magnitude larger than current binary pulsar bounds. 



1. Motivation 

One of the primary goals of the Laser Interferometer Space Antenna (LISA) is to study possible 
deviations from General Relativity (GR) in the strong field, where gravity is strong, speeds are 
large and binary pulsar or Solar System experiments are inapplicable. One route to such tests 
has been to study deviations from the Kerr metric within GR, through the determination of 
the multipolar structure of the gravitational field. In GR, however, the Kerr metric is known 
to be the only physically-reasonable representation of the exterior gravitational field of spinning 
compact objects, provided cosmic censorship and causality holds. 

Another route to test GR is to consider alternative theories of gravity and study the imprint 
these leave on gravitational waves (GWs) emitted by LISA sources. Here we report on the 
first consistent study of GW generation in a well-motivated alternative theory of gravity by the 



inspiral of a small compact object into a spinning, supermassive BHs (a so-called extreme-mass 
ratio inspiral or EMRI) [l]. 

Due to the plethora of alternative theories available, it becomes difficult to justify the choice 
of one theory over another. For this purpose, it is convenient to propose certain criteria for a 
theory to be a reasonable candidate to test GR with LISA [2j: 

• Metric theories of gravity: the gravitational field is described by a symmetric metric tensor 
that satisfies the weak-equivalence principle [3]. 

• Weak-field consistency: reduction to GR for weak gravitational fields and small velocities, 
such that experimental and observational tests are passed. 

• Strong-field inconsistency: deviations from GR in the dynamical strong-field, where gravity 
is strong and speeds are large. 

The weak-field consistency criterion also implies the existence and stability of physical solutions, 
such as the Newtonian limit of the Schwarzschild metric to describe physics in the Solar System. 
Additional criteria can also be proposed, such as well-posedness of the initial- value formulation, 
the existence of a well-defined and complete relativistic action, and some type of motivation 
from fundamental physics. These additional criteria, however, might be too stringent since only 
GR has been sufficiently studied to determine whether they are satisfied. 

A modification to GR that passes all the criteria in the itemized list above is Chern-Simons 
(CS) modified gravity [4J. In this 4D theory, the Einstein-Hilbert action is modified through 
the addition of the product of a dynamical scalar field and the Pontryagin density. Since the 
scalar field couples only to the quadratic curvature correction to the action, the weak-equivalence 
principle remains satisfied and the theory is indeed metric. Due to its quadratic curvature nature, 
this theory is also weak-field consistent but strong- field inconsistent. Stable solutions, like the 
Schwarzschild metric and a modified Kerr metric [5j, have also been found, although there are 
only qualitative arguments that suggest Cauchy well-posedness [6]. As for motivation, the CS 
modification arises generically and unavoidingly in the low-energy effective limit of string theory, 
both in the perturbative sector due to the Green-Schwarz anomaly-canceling mechanism [7J, 
and in the non-perturbative sector due to D-instanton charges and duality symmetries [8J. 
The CS correction can also be shown to arise in loop quantum gravity, when the Barbero- 
Immirzi parameter is promoted to a field in the presence of fermions in]- In fact, one can 
show generically (eg. in cosmological inffation) that the CS correction is one of a few terms 
that generically arise in effective field theories when one expands the action in the curvature 
tensor [TT] . 

Two different version of CS modified gravity exist, a non-dynamical and a dynamical one, 
where the former has been strongly constrained by theoretical and experimental considerations, 
while the latter has only bee weakly constrained by binary pulsar observations. The non- 
dynamical version assumes the scalar field is a priori prescribed and non-evolving, which 
forces what would have been an evolution equation for the scalar field into an additional 
constraint on the solution space of the theory. This additional constraint has been shown to 
be too restrictive, essentially disallowing spinning black hole (BH) solutions [I2] and forbidding 
perturbations of non-spinning BHs [13]. From an observational standpoint, the non-dynamical 
version has been strongly constrained in the neighbourhood of neutron stars through binary 
pulsar observations [14J; GWs would allow a constrain in regions of spacetime between source 
and observer due to corrections to GW propagation [151 [E] • 

Dynamical CS gravity has been only weakly constrained by binary pulsar observations [5j. 
Recently, |17 j studied perturbations of non-spinning, Schwarzschild BHs in the dynamical theory 
and found certain instabilities, which they used to argue a strong constraint could be placed 
due to the observationally confirmed existence of rotating BHs. This result, however, is strongly 
dependent on the assumption that there is no background scalar field to which the metric 



perturbations could couple to. Such an assumption is not valid for rotating BHs, even in the 
limit of small rotation, as found in [5]. Then, the work of [TT] can be considered as a good 
motivation to explore in more detail oscillations of BHs in the dynamical theory. As a first step, 
one would study the behaviour of linear oscillations for rotating BHs (which has not been done 
yet as it can be prohibitively difficult). Then, in the case that instabilities are still present, one 
could ask whether non-linear effects can suppress the linear instabilities. If the answer is in the 
negative, one can then try to put constraints on the theory based on astrophysical observations. 

We here study EMRIs in dynamical CS modified gravity. We first show that point-particles 
follow geodesies in this theory as in GR, although the supermassive BH background is modified. 
Such a modification to geodesic motion leads to corrections in the generation of GWs that lead 
to a dephasing between CS waves and those expected in GR. This dephasing will not prevent 
detection of GWs with LISA, but instead it will bias the estimation of parameters, leading to an 
uncontrolled systematic error. We conclude with a study of radiation-reaction in dynamical CS 
modified gravity and we proof that to leading order this is the same as in GR, except perhaps for 
subdominant energy-momentum emission by the scalar field. The inclusion of radiation-reaction 
will lead to stronger GW modifications that will break the degeneracy between the CS coupling 
parameter and the system parameters, thus allowing for a more stringent test of the dynamical 
theory relative to the binary pulsar constraint. 

The remainder of this Proceedings describes the above results in more detail: Section [2] 
presents the basics of the theory; Section [3] describes the generation of GWs; Section S] shows 
the difference between GR and CS GWs and discusses the inclusion of radiation-reaction; 
Section [5] concludes and points to future work. We employ the conventions in Misner, Thorne 
and Wheeler [3], described in more detail in [1]. 

2. Chern-Simons Basics 

We here summarize the basics of CS modified gravity and we refer the reader to the review 
paper [i\ for further details. 



2.1. Field Equations 

Consider then the Lagrangian density 



(1) 



where the gravitational constant is k = (IGvrG)"^; the coupling constant associated with the 
gravitational CS correction is a, while that associated with the strength of the CS scalar field 
'd is (3. The parity- violating Pontryagin density *RR : 



^^^^^^75^^/2, where the asterisk 



denotes the dual tensor, constructed using the antisymmetric Levi-Civita tensor e' 
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potential term V{'d) will be assumed to vanish, as it does in string theory due to shift symmetry. 
The last term in ([T]) represents the Lagrangian density for external matter degrees of freedom. 
Variation of the action with respect to all degrees of freedom leads to the CS field equations: 
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where T^"" is a matter stress-energy tensor, T^^y' is the stress-energy of the CS scalar field and 
C^^ is the so-called C-tensor, both of which are given by 

1 



(V'^??)(V,i?) 



As we are here interested in the dynamical version of dynamical CS modified gravity, we 
restrict attention to coupling constants /? / 0. If /? = 0, one would recover non-dynamical CS 
theory and the equation of motion for ?? would reduce to the Pontryagin constraint *RR = 0, 
which over-constrains the field equations [13]. We consider the dynamical theory, where we 
define ^ := o? /{k[3) as an expansion parameter for the so-called small-coupling approximation 
(see [5] for details), in which a modified Kerr solution is obtained. 

2.2. Gravitational Waves in CS theory 

Let us decompose the metric tensor as g^^ = g^j^ + e h^^ -l-C'(e^), where g^^ is some background 
metric, /i^^ is a metric perturbation, which in the radiation-zone (many GW wavelengths away 
from the center of mass) can be associated with GWs, and e is a book-keeping perturbation 
parameter associated with the GW strength. 

The CS field equations can now be expanded in e. In the radiation zone, however, we can 
further simplify the analysis by focusing on a flat background, g^^ = rj^^. We then find 

AC 

am = -|e2 e-'3M-/,^,,^^/,;7,<5]^ + 0(6^) , (6) 

where T^^ = T^y — {l/2)7]^j^T is the trace-reversed version of the total stress-energy tensor 
T^^ = TJ^j^^+Tj^\ = rj''''d^d,y is the flat-space D'Alembertian and e'"'°''^ is the flat-space Levi- 
Civita tensor density (not to be confused with e). We work in Lorenz gauge, (/i'^'^ — ^hg^^'^^ |^ = 0, 
where h = g^'^hfj,^, is the trace of /i^^, with respect to g^^, and we restrict the gauge further to a 
transverse-traceless one hy h = 0. 

The GW generation formula are then identical to the GR ones to leading order in e, as one 
can see in ([5]). This is because the CS correction is proportional both to the scalar field and 
the metric perturbation, where the former is sourced by the Pontryagin density [see ([6])], which 
itself is of 0{e'^) in the radiation zone. One can make this argument more formally accurate by 
decomposing the scalar field itself into a background plus a perturbation and then performing a 
multiple-scale perturbative expansion, but we relegate such details to \X\. Furthermore, one can 
show that there are only two independent GW polarizations in dynamical CS gravity, as is the 
case in GR, in spite of the extra scalar field present in the theory. This can be done by studying 
the geodesic deviation equation or plane wave solutions, as found in [UIB]. 

3. EMRIs in CS Gravity 

We now consider the generation of GWs by EMRIs in dynamical CS gravity. We shall 
approximate these waves in the so-called semi-relativistic approximation |18j , where the motion 
is assumed geodesic and GWs are assumed to propagate in flat spacetime. This simplifying 
assumptions have recently been shown to be accurate relative to more precise black-hole 
perturbation theory calculations [19]. We begin with a description of the background upon 
which geodesies are evolved. We then show that point-particle orbits can truly be described with 
geodesies and explicitly derive the latter. We end by presenting the GW generation formula in 
the semi-relativistic approximation. 

3.1. Background Geometry 

Via the small-coupling approximation, complimented by a slow-rotation approximation, one can 
solve for the exterior gravitational field of a rotating BH in dynamical CS modified gravity. 



The line element is given by ds'^ = ds|,,, + 5^a/(4r^) [l + 12M/(7r) + 27M^/{Wr^)] sin^ edtd(p, 
where dsxerr ^^^^ element for the Kerr metric and we have used Boyer-Lindquist type 

coordinates [5]. 

The CS modified metric remains stationary and axisymmetric, allowing for the definition 
of multipole moments appropriately. The multipolar structure of the modified metric remains 
completely determined by only two moments (no-hair or two-hair theorem): the mass monopole 
and the current dipole. The relation, however, between these two moments and higher order 
ones is modified from the GR expectation at multipole £ > 4. In spite of the CS scalar field, 
the dynamical theory satisfies the no-hair theorem, as the scalar field is fully determined by the 
geometry. 



3.2. Motion of Massive Bodies 

One can show that test-bodies follow geodesies in CS modified gravity [Ij. This is because the 
divergence of the field equations ([2]) leads to ^ for the non-matter degrees of freedom, forcing 
the divergence of the matter stress energy tensor to vanish. This then proves dynamical CS 
modified gravity satisfies the weak-equivalence principle, thus rendering it a metric theory as we 
anticipated in the Introduction. 

The motion of small compact object with an extreme-mass ratio can then be modeled, 
to zeroth order, by geodesies in the modified Kerr background of the previous section. As 
the background remains stationary and axisymmetric, there still exist timelike and azimuthal, 
commuting Killing vectors that define a conserved energy E and angular momentum L per 
unit mass. Moreover, the modified background also possesses an additional constant of the 
motion (the Carter constant, Q) associated with a Killing tensor. Using these constants 
and the orthonormality condition for timelike geodesies, we can derive the geodesic equations 

= x^gj-j, -|- Sxcs, where an overhead dot stands for partial time derivative and x'^ = [t, r, 9, (p] 
are Boyer-Lindquist coordinates (where the polar angle 9 is not to be confused with the CS scalar 
field i}). The correction factor fe^g = [L,2EL,0, -E]5g'^'' , where dg^"" = ^a/{ll2r^ f)i70r^ + 
UOr M + 189M2) and / := 1 - 2M/r. 

The geodesic equations allow one to compute some physical observables that waveforms are 
sensitive to. One such observable is the innermost-stable circular orbit (ISCO) location, which 
is CS shifted by [5]: R^gco = 6M =F 4\/6a/3 - 7a^/{l8M) ± 77\/6a^/(5184M^), where the 
upper (lower) signs correspond to co-rotating (counter-rotating) geodesies. Notice that the CS 
correction works against the spin effects. Other observables of interest are the fundamental 
orbital frequencies associated with orbital motion [20l[2l]. In CS theory, the temporal, radial 
and azimuthal frequencies are modified in an unilluminating way that we avoid here [T] , if one 
employs the time coordinate introduced in [22j. Therefore, CS geodesies correspond to GR 
geodesies with different fundamental frequencies, which maps to different system parameters. 



3.3. GWsforEMRIs 

The first simplification of the semi-relativistic approximation (that of modeling orbital motion 
as geodesic) has already been discussed in the previous section, and so we here concentrate on 
the second semi-relativistic simplification: GWs propagating on fiat spacetime. We implement 
this approximation by modeling the solution of the linearized GW generation formula in [5] 
via a multipolar expansion (see eg. [23]). In particulat, we employ a quadrupole-octopole 
approximation, where the plus- and cross-polarizations of the waveforms are given by 



1 i 



lij 2n^Sfiij -\- I 
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where r is the flat-space distance from the source to the observer. 



[l,x7r] 



are 



polarization tensors, and lij, Sijk and lijk are the symmetric and trace- free projections of the 



mass quadrupole, current octopole and mass octopole respectively, defined as integrals of the 
product of point-particle stress-energy tensor with the point-particle trajectories (see eg. [23] for 
explicit expressions). 

An important detail in the above implementation is the relation between the coordinates used 
in the geodesic evolution (Boyer-Lindquist) and those employed in the metric reconstruction 
(Cartesian). We here choose to identify Boyer-Lindquist coordinates with flat-space spherical 
coordinates, such that {x,y,z) = r(sin cos sin sin 0, cos 6*). Another important detail is 
that the above approximation is truly a slow-motion one, which implies waveforms are accurate 
provided the pericenter orbital velocity is much smaller than the speed of light. In spite of 
this, the semi-relativistic approximation does capture the correct qualitative behaviour of EMRI 
orbits and waveforms. 

4. Testing CS theory with LISA and Radiation-Reaction 

We here describe some of the differences between a GR and a CS EMRI waveform. 
For this, we shall concentrate on the following test orbits: (A) {cL,(,,rp^^i,e,9i^J = 
(0.1M,0.1M^12M,0.2,0.1), (B) {a,C,r^^^„e,e,J = (0.2 M, 0.2M^ 8 M, 0.4, 0.2) , (C) 
(^) ''pori) ^1 ^inc) = (0.4M, 0.4M^, 6 M, 0.6, 0.3) , where r^^^^ stands for pericenter distance, e for 
orbital eccentricity and 9^^^ is related to the inclination angle. Waveforms are measured by z-axis 
observers at a distance of 8 kpc and test orbits have been ordered from least to most relativistic. 
All orbits assume the central BH has M = M,, where M, = 4.5 x 10^ Mq is approximately the 
mass of the BH at the center of the Milky Way, while the SCO has m = 35Mq, which leads to 
a mass ratio fj, ~ 7.8 x 10~®. We conclude with a discussion of radiation-reaction in CS theory. 

4.I. Trajectory and Waveform Dephasing 

The trajectories are obtained by solving the geodesic equations discussed above for a total time 
of r = 5 X lO^M, obtaining on the order of 10^-10^ cycles, (see eg. [T] for details). Geodesies 
are initialized with the same orbital data when performing GR and CS evolutions. That is, 
given some set of conserved quantities {E,L,Q), one can derive three initial orbital parameters 
{p,e,9.^J, but these will differ in GR and CS due to the latter's geodesic corrections. We here 
consider orbits with the same orbital parameters, instead of the same constants of motion. 

The orbits chosen present rather generic behaviour with a stage of zoom-whirl, where the 
particle whirls violently for several cycles close to pericenter and then zooms out to large radius 
only to return again close to pericenter, during all of which there is both in-plane and out-of- 
the-plane precession. The left panel of [T] shows the projection of the orbital trajectory onto 
the x-y plane (orbit C) assuming a CS modified Kerr background BH (black line) and a Kerr 
background (light gray line). We present only the last 17500M of geodesic evolution, during 
which we can clearly observe that the trajectories have dephased significantly . We find that 
orbit C presents the biggest dephasing, followed by B and then A. 

The right-panel of [U shows the difference in the plus-polarized waveforms in CS gravity and 
GR. As a reference, the maximum magnitude of the GR GW polarization is approximately 
l^+l < 8 X 10"^^ (orbit A), \h^\ < 1.25 x lO'^^ (orbit B), and < 2.25 x lO"^*^ (orbit 
C). Observe that by the end of the simulation (~ 128 days of data using M = M,), orbit A 
waveforms have dephased by 0.3% , while orbit B waveforms have dephased by 16%, and orbit C 
waveforms by 90% relative to the maximum amplitude of the respective GR waveforms. Similar 
behaviour is observed for the other polarization. 

4-2. Radiation- Reaction in Dynamical CS Gravity 

Radiation-reaction can be incorporated via a radiative, adiabatic approximation |22l [25l [26], 
where the constants of the motion are evolved using balance laws that relate orbital energy loss 
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Figure 1. Left: Two-dimensional projection onto the x-y plane of orbit C in a Kerr background 
(light gray) and in the CS corrected Kerr metric (black). Observe that a large dephasing has 
built up by the end of the evolution. Right: Plus-polarized GW dephasing as a function of 
time. The top panel shows results for orbit A, the middle one for B and the bottom one for C. 
Observe that after 2/5 of the evolution, the waveform has dephased significantly. 



to GW energy emission. Such an approximation is valid provided the change in the constant of 
the motion occurs on a much larger timescale relative to the orbital one. 

The balance laws can be obtained by studying the effective stress-energy tensor associated 
with GWs in the framework of the short-wavelength approximation |27[|28j . We thus decompose 
both the geometry and the CS scalar field into a slowly- or non-oscillating background and 
a rapidly-oscillating perturbation. The linear part of the modified field equations describe 
the evolution of the perturbation, while the averaged second-order part describes how the 
background is modified by GWs. This later part serves as an effective GW stress-energy tensor, 
T^^ (the Isaacson tensor), 

+^ (< ^'^c^A^M > + < '-^'c^A^M >)} (8) 

where the superscript preceding any given quantity stands for the perturbative order of that 
quantity with respect to metric perturbations. Restricting attention to the transverse-traceless 
gauge and using properties of the averaging procedure, one can show that all terms associated 
with the CS correction in [H] identically vanish, yielding the same expressions for the Isaacson 
tensor as in GR. 

The above result proofs that orbital backreaction in dynamical CS modified gravity is similar 
but not identical to that expected in GR. First, in addition to GW emission, orbital energy is 
also lost via scalar field emission, although the latter is a subleading effect. Second, even if one 
neglects such scalar field emission, GW emission depends on time-derivatives of the waveforms, 
which we showed in the previous subsection are different in CS gravity. 

4-3. A possible LISA Implementation of a CS Test 

What is the effect of the CS correction in LISA data analysis? As shown through the semi- 
relativistic waveforms, GR and CS waveforms dephase given sufficient time, but such dephasing 
will not prevent detection. Instead such dephasing will only affect parameter estimation. That 



is, the CS correction modifies the fundamental geodesic frequencies, and thus, a CS geodesic is 
equivalent to a GR geodesic with different orbital parameters. 

How big should the CS coupling parameter be to induce a sufficiently large dephasing 
that would contaminate parameter estimation? [2] shows the average of a normalized dephasing 
measure as a function of time for orbit C with different ^ (see caption) . This dephasing measure 
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Figure 2. Normalized dephasing measure ©(gr, cs) as a function of time for orbit C with: 
^ = (solid black), i = 0.2M^ (dotted blue), ^ = O.IM^ (dashed red), ^ = O.OBM^ (dot- 

dashed green), and ^ = 0.025M^ (dot-dot-dashed violet). Observe that a significant dephasing 
is found even for ^ = 10^^. 



is nothing but the integrand of the overlap with averaged beam-pattern functions in the time- 
domain: O(gR)Cs) := [h'^^h'^ + /ix^/ix^] normalized by OAr(GR,cs) := y^O(GR, gr) 0(cs,cs). 
This figures shows that the CS correction could lead to a significant dephasing, and thus bias in 
parameter estimation, even for a CS parameter of 0(10"^) after only 4 months of data. 

The inclusion of radiation-reaction is crucial for breaking the degeneracy between the CS 
correction ^ and the system parameters, i.e. to distinguish between a GR and a CS geodesic. 
One expects that radiation-reaction will not only break this degeneracy, but also lead to larger 
deviation from the GR expectation. A rough estimate of the accuracy to which CS gravity 
could be constrained via a LISA observation is ^^^^ < 10^ km (5/10^®)"'^^^ (M/M,), where S is 
the accuracy to which ^ can be measured, which depends on the integration time, the signal- 
to-noise ratio, the type of orbit considered and how much radiation-reaction affects the orbit. 
Notice that intermediate-mass ratio inspirals (IMRIs) are favored over EMRIs. This result is to 
be compared with the binary pulsar constrained ^^/'^ < 10^ km [5]. We then see that an IMRI 
with M = IO^Mq could pi constraint two-orders of magnitude more stringent than the 
binary pulsar one. Moreover, a GW test can constrain the dynamical behaviour of the theory 
in the neighbourhood of BHs, which is simply not possible with binary pulsar observations. 

5. Conclusions 

Dynamical CS gravity is a viable theory that can only be constrained dynamically via GW 
observations. We have constructed geodesies in the modified theory and derived radiation- 
reaction formula that should allow for the construction of CS-modified EMRI waveforms. Such 
waveforms could be used in LISA data analysis pipelines to test the theory at least two-orders 
of magnitude better than current binary pulsar observations. 
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